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Abstract 



We present exact calculations of reliability polynomials R(G,p) for lattice strips G of 
fixed widths L y < 4 and arbitrarily great length L x with various boundary conditions. 
We introduce the notion of a reliability per vertex, r({G},p) = lim^y^^ ^(G^p) 1 /'^' 
where |V| denotes the number of vertices in G and {G} denotes the formal limit 
limiy[_ K)0 G. We calculate this exactly for various families of graphs. We also study 
the zeros of R(G,p) in the complex p plane and determine exactly the asymptotic 
accumulation set of these zeros B, across which r({G}) is nonanalytic. 
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1 Introduction 



With the development of progressively larger communications networks such as those in- 
volving telecommunications and those connecting computers, from large local area networks 
to the internet, the analysis of the reliability of these networks has become an increasingly 
important area of study. Because many types of failures of communication links are random, 
this study falls within the area of statistical physics, as well as engineering and, as will be 
evident, mathematical graph theory. A common approach to this study is to factorize the 
reliability into a factor describing the components such as message routers and computers at 
various nodes, on the one hand, and the links between these nodes, on the other hand. One 
often assumes that nodal network elements like message routers and computers have a certain 
probability p no de of (normal) operation. Since failures of these nodal components generally 
occur independently, the factor in the overall reliability due to these elements is then p" ode , 
where n is the number of nodes in the network. The other factor contributing to the overall 
network reliability involves the connecting links on the network and hence the connectivity 
structure of the network. As mentioned above, many failures in these connecting links are 
also random and independent and may thus be treated in a probabilistic manner in which 
one describes a given link as having a probability p of operating and thus a probability \ —p 
of failing. (Examples of nonrandom failures include earthquakes that affect only one portion 
of a network; these are not considered here.) The network is thus represented as a connected 
graph G = (V, E) with vertex (= node) set V and edge (= bond) set E, in which each 
vertex represents either a message router or switch or an endpoint such as a computer or 
terminal (or telephone in the telecommunications example) and each edge represents the link 
between the nodes. We shall be interested here in a commonly used simplified model of the 
network in which each nodal element and link have, respectively, the same fixed probabilities 
Pnode and p of operation with p no d e and p lying in the interval [0, 1]. We shall denote the 
number of vertices and edges in G as \ V\ = n and \E\. In the graphical representation of the 
network, one regards a normally operating communications link as an edge that is present 
(with probability p), and a malfunctioning communications link as an edge that is absent 
(with probability 1 — p) , the presence or absence being random and uncorrelated. Defining a 
successfully operating network as one in which all of links and nodes are operating normally, 
it follows that the node reliability is, as noted above, p^ de , independent of the connectivity of 
the network. The factor due to the operation of the links is much more difficult to calculate. 
It is therefore customary, in modelling the reliability behavior of networks, to separate out 
the node-reliability factor and concentrate on the contribution due to the structure of the 
network. 

The all-terminal reliability polynomial R(G,p) is defined as the probability that there 
is an operating communications link between any two nodes in the network, i.e. any two 
vertices in the set V are connected (by a path consisting of edges that are present) [l]-[5]. 
In passing, we note that one can also study /c-terminal reliability probability polynomials, 
but we shall restrict ourselves here to the all-terminal reliability polynomial and hence shall 
henceforth omit the qualifier "all-terminal" in the notation. The contributions to R(G,p) 
arise from the sum of connected spanning subgraphs of G. Here a spanning subgraph is a 
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subgraph G' = (V, E') with E' C E, i.e. a graph with the same vertex set and a subset of the 
edge set of G. Denote a connected spanning subgraph as H = (V,E H ). Each contribution 
to R(G,p) is weighted by a factor p\ E "\ reflecting the probability that the \E H \ edges are 
present, multiplied by a factor (1 — p)\ E \~\ E n\ reflecting the probability that the edges in the 
complement set {e G E; e ^ Eh}, of cardinality \E\ — \Eh\, are absent. That is, 

R(G,p) = J2 P ]Eh \ 1 -p) IEHEh1 . (1.1) 

HCG 

Some pioneering studies of network reliability include [l]-[4]; a review is [5]. The calculation 
of R(G ) p) for an arbitrary (connected) graph G has been proved to be #-P complete [5, 6]; 
that is, roughly speaking, for a generic network G, the time required to calculate R(G,p) 
grows exponentially with the size of the network. In the face of this difficulty, network 
designers have often relied upon upper and lower bounds on R(G,p) [7]- [13]. It is clearly of 
interest if one can obtain exact solutions for R{G,p) for some families of graphs G. 

In this paper, we shall present a number of such exact solutions for reliability polyno- 
mials, for recursive families of graphs G. We have previously presented such solutions for a 
specific recursive family of graphs in [14]. A recursive family of graphs is one in which one 
constructs successive members of the family in a recursive manner starting from an initial 
member. Recursive families of graphs that are of particular interest here are strips of regu- 
lar lattices of a given width L y vertices and arbitrarily great length L x vertices with some 
specified transverse and longitudinal boundary conditions. We shall envision these strips 
as extending in the horizontal (= longitudinal, x) direction and having transverse extent in 
the vertical, y, direction. To see that these form recursive families, one can picture a strip, 
say of the square lattice, of length L x and width L y as being formed by gluing on a column 
of squares of height L y to the strip of length L x — 1. The boundary conditions that we 
shall consider include free (FBC), periodic (PBC), and twisted periodic (TPBC), by which 
is meant that the longitudinal ends are identified with a reversal of orientation. We shall 
denote the various combinations as (i) (EBC y , FBC X ) = free, (ii) (FBC y , PBC X ) = cyclic, 
(iii) (FBC y ,TPBC x ) = Mdbius, (iv) (PBC y ,FBC x ) = cylindrical, (v) (PBC y ,PBC x ) = 
torus, and (vi) (PBC y , TPBC X ) = Klein bottle. For an arbitrary graph G, one defines the 
degree di of a vertex Vi as the number of edges connected to it. One then defines the maxi- 
mum degree as A = max v . e y di. A A-regular graph is one in which all of the vertices have 
the same degree A. Some strip graphs, such as those with toroidal or Klein bottle boundary 
conditions, are A-regular graphs. For strip graphs that are not A-regular, it will still be 
useful to define an average or effective vertex degree in the limit of infinite length, namely 

Secondly, we shall introduce a definition of the reliability per node in the limit of many 
nodes, 

r({G},p)= lim RiCp) 1 ^ (1.3) 
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where {G} denotes the formal limit of a given family, hin^^oo G. We calculate this function 
exactly for a number of families and study its properties. As we shall show, this function 
provides a convenient quantitative measure of the reliability properties of large networks. 

A third part of our work involves the calculation of the zeros of R(G,p) in the complex 
p plane and the determination of their asymptotic continuous accumulation set, denoted B 
as the number of vertices \V\ — > oo (in practice, for our strips, this is equivalent to the limit 
of infinite length, L x — > oo). As we shall show, the function r({G}) is nonanalytic across 
the locus B. We define p c as the maximal point where B intersects (or coincides with) the 
real axis. For the infinite-length limits of some lattice strips, the locus B does not cross the 
real axis, so no p c is defined. In certain of these cases, there are complex-conjugate arcs on 
B whose endpoints at p en d, p* end are very close to the real axis; in these cases, it is useful to 
define a (p c )eff — Re(Pend)- Although the width L y — > oo limit of loci B for the infinite- 
length lattice strips is not necessarily the same as the continuous accumulation locus of the 
zeros of the reliability polynomial for the two-dimensional thermodynamic limit L x — > oo, 
L y — > oo with L y /L x fixed to a finite nonzero constant, our exact results on infinite-length 
limits of lattice strip graphs can give some insight into plausible behavior of the locus B for 
the thermodynamic limits of the corresponding 2D lattices. We shall discuss this further 
below. 

There are several motivations for our study. Clearly, new exact calculations of reliability 
polynomials are of interest from the point of view of both the statistical physics of network 
theory and of mathematical graph theory. The generalization of p from the interval p G [0, 1] 
to p G C is clearly necessary in order to analyze the zeros of the reliability polynomial and 
their accumulation set as \ V\ — > oo. Our results also further illuminate the fascinating variety 
of applications of the Tutte polynomial, or equivalently, the Potts model partition function. 
Our focus here is thus somewhat more abstract than in a specific engineering context, where 
one is trying to design the most reliable network subject to various constraints, such as cost. 



2 Some General Properties of the Reliability Polyno- 
mial 

In this section we review some general properties of the reliability polynomial that will be 
used in the rest of the paper. As is evident from its definition, the reliability polynomial has 
the properties that 

R(G,0) = 0, R(G, 1) = 1, (2.1) 
IfpG[0,l], then R(G,p) G [0, 1] (2.2) 

and 

^^>0 for pG[0,l]. (2.3) 
dp 

In general, an edge can be classified as (i) a loop, i.e., an edge that connects a vertex 
to itself, (ii) a bridge (= co-loop) with the property that if it is deleted, this increases the 
number of connected components of the graph by one (so that, if G is connected, deleting 
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the bridge breaks G into two disjoint parts), or (iii) neither a loop nor a bridge. We shall use 
the standard mathematical notation that G/e means the graph obtained by deleting e and 
identifying the two vertices that were connected by e (called contracting G on e) and G — e 
means the graph obtained by deleting e. Now if e is a loop, e = ee, then the probability that 
all terminals are connected with each other is independent of ee, so 

R(G,p) = R(G/e e ,p) . (2.4) 

If e = e& is a bridge, then the probability that it is present is an overall multiplicative factor 
in R, so 

R(G,p)=pR(G/e b ,p) . (2.5) 
If e is neither a loop nor a bridge, then R(G,p) satisfies the deletion-contraction relation 

R(G,p) = (l-p)R(G-e,p)+pR(G/e,p) . (2.6) 

That is, if the edge e is not present, then R(G,p) has the value given by the reliability 
polynomial of a graph G — e times the factor (1 — p); if the edge is present, then R(G,p) has 
the same value as the reliability polynomial of a graph contracted on e, multiplied by p, and 
since these two possibilities (e not present or present) are mutually exclusive and exhaust 
the possibilities, the full polynomial R(G,p) is the sum of these two possibilities. 

We next comment on some general properties of our r({G},p) function. We recall that 
a tree graph T n is a connected graph with n vertices and no circuits. Since R(T n ,p) = p n_1 , 
it follows that r({T},p) = p, where {T} = lim^oo T n . However, in general, r differs from p 
for network topologies that are more complicated than that of a tree graph. The r function 
can give insight into the reliability behavior of large networks. From obvious arguments, 
it follows that any reliability function, whether the /c-terminal reliability function Rk(G,p) 
or the all-terminal reliability Ra{G,p), is a monotonically increasing function of p on the 
relevant interval p G [0, 1]. It follows that r({G},p) is also a monotonically nondecreasing 
function of p on this interval. 

From the basic definition (1.1) one can generalize p from the physical interval < p < 1 
to arbitrary real or, indeed, complex p. This generalization is necessary when one calculates 
the zeros of R(G,p) in the complex p plane. Using our calculations of R(G,p) for a variety 
of families of lattice strip graphs, we shall determine exactly the continuous accumulation 
set of the zeros of R(G,p) in the complex p plane as \V\ — > oo for each family of graphs 
G. We shall denote this accumulation set as B. In this context we note that in 1992 Brown 
and Colbourn conjectured that for an arbitrary connected graph G, all of the zeros of the 
reliability polynomial R(G,p) lie in the disk \p— 1| < 1 [15] (some related work is [16, 17, 18]). 
Recently, Royle and Sokal have reported counterexamples that show the Brown-Colbourn 
conjecture to be false [19], although the maximal values of \p — 1| found so far are only 
slightly greater than unity. Extending the counterexamples found in [19], we have also found 
several recursive families of graphs for which the conjecture is false. 

One can also address a related question for the continuous asymptotic accumulation set 
B: Consider a recursive family of arbitrary connected graphs G = G(V, E) and, as above, 
let B be the asymptotic accumulation set of zeros of R(G,p) in the limit |V| — > oo. Is it 
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true that all points on B are contained in the disk disk \p — 1| < 1? Clearly, although the 
Brown-Colbourn conjecture is false, the answer to the above question could still be yes since 
the zeros of R(G,p) that lie outside the disk \p — 1| might, a priori, not accumulate to form 
part of B. However, we shall show that the answer to this question is no, i.e. we shall exhibit 
recursive families of graphs for which B extends outside the disk \p — 1|. Using our exact 
results, we shall address a number of other interesting questions concerning B. It is also of 
interest to note that Brown and Colbourn proved that all of the real roots of the reliability 
polynomial lie in the set {0} U (1,2] [15]. 

Since we shall present calculations for strip graphs with free and periodic or twisted 
periodic longitudinal boundary conditions, we observe the following inequality, where BC y 
refers to any given transverse boundary condition 

Theorem 2.1 For a strip graph G s of length L x and width L y vertices, 

R(G s ,L y x L x ,BC y ,(T)PBCx,p) > R(G S , L y x L x , BC y , FBC x ,p) (2.7) 

Proof These inequalities follow because for a given type of lattice strip graph G s with 
a given length L x and width L y , the strip with periodic or twisted periodic longitudinal 
boundary conditions has a greater connectivity than strip with free longitudinal boundary 
conditions. In particular, for a given G s , if all of the transverse edges at a fixed longitudinal 
location are absent, this does not produce two disjoint graphs for a cyclic or Mobius strip 
but does for the free strip. □. 



3 Relations Between the Reliability Polynomial and 
the Tutte Polynomial and Potts Model Partition Func- 
tion 

We recall that the Tutte polynomial of a graph G = (V, E) is [20]- [27] 

T{G^y)=Y,^-^ kiG ' ) - k{G \y-^ iG ' ) (3-1) 

G'CG 

where k(G') and c(G') denote the number of components and linearly independent circuits of 
G', with c(G') = \E'\+k(G') — \V\. Since we deal only with connected graphs here, k(G) = 1. 
(We follow standard notation for the variables x and y of the Tutte polynomial and caution 
the reader that these have no connection with the x (longitudinal) and y (transverse) axes of 
the strip graphs.) The reliability polynomial R(G,p) can be expressed in terms of the Tutte 
polynomial T(G, x, y) as [5] 

R(G,p)=p\ v \-\l- P r\-\ v ^T(G,l,- [ ^) . (3.2) 
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We shall use this relation to calculate reliability polynomials for various families of graphs 
from our previous calculations of the Tutte polynomials for these graphs [28]- [38]. We have 
also used an iterative application of the deletion-contraction relation (2.6) to calculate reli- 
ability polynomials directly. 

The Tutte polynomial is equivalent to the partition function of the Potts spin model. On 
a lattice, or more generally, on a graph G, at temperature T, this model is defined by the 
partition function [39, 40] 

Z(G,q,v) = J2 e ~ m ( 3 - 3 ) 

Wi} 

with the Hamiltonian 

n = -jJ2 6 ^ ( 3 - 4 ) 

where a, = 1, ...,q are the spin variables on each vertex i e G; f3 = (A; B T) _1 ; (ij) denotes 
pairs of adjacent vertices; and we use the notation 

K = (3J , a = e K , v = a - 1 . (3.5) 

This partition function can be written as [41] 

Z(G,q,v)= g fc(G V (G0 • (3-6) 

G'CG 

This formula enables one to generalize q from Z + to R or, indeed, C. From it one also 
directly deduces the equivalence 

Z(G, q, v) = (x- l) k ^ G \y - l)l y lT(G, x, y) (3.7) 

where 

x = 1 + - (3.8) 

v 

and 

y = v + 1 = e K (3.9) 

so that 

q = (x-l)(y-l) . (3.10) 

Combining (3.2) and (3.7), we have, for a connected graph G, the relation between the 
reliability polynomial and the Potts model partition function 

R(G, V ) = (l- V f \ lini i q -i Z (G,q,v = -2-) ■ (3.11) 
<?^o i — p 

Note that Z(G,q,v) always has an overall factor of q, which cancels the factor of q~ x in 
(3.11). The valuation y — 1/(1 — p) in (3.2), or equivalently, v — p/(l — p) in (3.11) can also 
be expressed as 

p = l- e - K = ^^. (3.12) 

Thus the physical range of p G [0, 1] for the network corresponds to the physical range of 
temperature for the q = Potts ferromagnet, with T = p — 1 and T = oo p — 0. 
The physical range of temperature for the Potts antiferromagnet would correspond to the 
unphysical interval — oo < p < for the network. 
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4 Structural Results 



4.1 General Recursive Families of Graphs 

A general form for the Tutte polynomial for the strip graphs considered here, or more 
generally, for recursively defined families of graphs G m comprised of m repeated subgraph 
units, is [28] 

T(G m ,q,v) = — — <*A*Gj) m ( 4 - L1 ) 

X 1 

where the terms X G ,j, the coefficients cc,j, and the total number Nz,g,\ depend on G through 
the type of lattice, its width, L y , and the boundary conditions, but not on the length. 
We observe the following general result: 

Theorem 4.1 Consider a recursive family of graphs G m comprised ofm repeated subgraph 
units. Then the reliability polynomial R(G m ,p) has the general form 

R(G m ,p) = J2 <-'<,;. J (li<;,r (4.1.2) 

where the coefficients c^gj an d terms fJ-o/s depend on G through the type of lattice, its 
width, L y , and the boundary conditions, but not on the length. 

Proof This is proved by combining (3.2) and (4.1.1). Clearly, the A*gj' s an d cr^/s 
can be determined from the A Gj 's and c G /s. □ The reduction of the form (4.1.1) to (4.1.2) is 
analogous to the reduction that we had found for chromatic polynomials of recursive families 
of graphs [42] -[50]. 

It is convenient to factor out a power of p and write (4.1.2) as 

R(G m ,p)=p a ^ m+a -° c W«oj) m (4- 1 - 3 ) 

j'=i 

with 

a G j=p- ai - a li G j (4-1.4) 

where again the a^a, i — 0,1, depend on G through the type of lattice, its width, L y , and 
the boundary conditions, but not on the length. Obviously, N^c,a = N^g,^- We find that 
for p = 1, for all of the families that we have considered, all except one of the N^c,a «gj's 
vanish, and the nonzero term, which can be taken to be the first, has the value unity: 

a G ,i = l, a G:j = 0, 2<j<N R:G , a for p = 1 (4.1.5) 

For the strips of the square and triangular lattices considered here, 

a hG = L y , G of sq,t type (4.1.6) 
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and for the strips of the honeycomb lattice, 

a ltG = 2Ly, G of he type (4.1.7) 

Note that for all of the strips of the square and triangular lattice with periodic or twisted 
periodic longitudinal boundary conditions, m = L x . For the corresponding strips with free 
longitudinal boundary conditions, it is convenient to set m — L x — 1; with this definition, it 
follows that for strips with both free and periodic (or twisted periodic) longitudinal boundary 
conditions, m is equal to length measured in units of edges, with L x being the length in terms 
of vertices. The situation with the honeycomb strips is different. For example, for a cyclic 
strip of the honeycomb lattice, envisioned as a brick lattice, L x = 2m; that is, if the strip is 
m bricks long, then the length, measured in terms of vertices, is 2m. 

Following our earlier nomenclature [42], we denote a term acj in (4.1.3) as leading (= 
dominant) if it has a magnitude greater than or equal to the magnitude of other envy's. 
In the limit n — > oo the leading aaj in R(G,p) determines the function r({G},p). The 
continuous locus B where r({G},p) is nonanalytic thus occurs where there is a switching 
of dominant aajs in R, respectively, and is the solution of the equation of degeneracy in 
magnitude 

\<xgj\ = Wg,A (4.1.8) 

between two dominant cc's. 



4.2 Cyclic Strips of Regular Lattices 

Here we will derive some structural properties of reliability polynomials for cyclic strips of 
regular lattices from corresponding properties of Tutte polynomials, using eq. (3.2). In [31] 
it was shown that for cyclic and Mobius strips of the square lattice of fixed width L y and 
arbitrary length L x (and also for cyclic strips of the triangular lattice) the coefficients Cj in 
the Tutte polynomial are polynomials in q with the property that for each degree d there is 
a unique polynomial, denoted S d \ Further, this was shown to be 

c«> = tW /2 /2) = Y,{-l) j ( 2d ~ j \ d - j (4.2.1) 
j=o \ 3 / 

where U n (x) is the Chebyshev polynomial of the second kind. A number of properties of 
these coefficients were derived in [31]. The first few of these coefficients are 

c (o) = x ; c (i) = q _ 1 ^ c (2) = q 2 _ 3q + 1 ; (42>2) 

c (3) = g 3 _ 5g 2 + Qq _ l _ ^.S) 

Following our earlier work [36], we define 

K {d) = c {d) + c (d-i) 
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= £(-i)''( 2d • % d ~ 3 ( 4 - 2 - 4 ) 

3=0 

and 

= g-V«0 = ^(_i)jY 2d - 1 " j \ d-i-i . (4.2.5) 
j'=o \ 3 J 

An important property of is that it always has the factor q; this is evident from (4.2.4). 
The first few of these polynomials are 

/t (1) = g, /t {2) = g(g - 2) , /t {3) = q(q - l)(q - 3) . (4.2.6) 

The following special results for q — 0, i.e., x = 1, will be useful: 

UmR {d) = (-l) d+1 d (4.2.7) 

c (d) = (_^d for 5 = o ( 4>2>8 ) 

and 

^ = ( _ 1} «w+i)rf for9 = (42 . 9) 

2 

so that, using (3.10) to express as a function of x and y, 

^ = (-!)"(„ -1)<1+W for x=l. (4.2.10) 

The terms \T,G s ,L y ,j that occur in (4.1.1) can be classified into sets, with the \T,G s ,L y ,j(x, y) 
in the cPth set being defined as those terms with coefficient S d \ In Ref. [31] the numbers 
of such terms, denoted n,T{L y , d), were calculated. Labelling the eigenvalues with coefficient 
as \T,G a ,L y ,d,j with 1 < j < n,T(L y ,d), the Tutte polynomial for a strip graph of length 
L x = m of a regular lattice of type G s can be written in the form 

L y n T (L y ,d) 

T[G s (L y xm;FBC y ,PBC y ),x,y] = -—J2 c(d) E (^) m ■ ( 4 - 2 - U ) 

The property that T(G, x, y) is a polynomial in x, which is evident from its definition (3.1), is 
not manifest in (4.2.11) because of the l/(x — 1) prefactor; however, this prefactor is always 
cancelled in the evaluation of the sum on the right-hand side of (4.2.11) for any particular 
type of lattice strip graph G m . Nevertheless, the presence of this prefactor means that if 
one calculates the reliability polynomial for a cyclic strip graph by evaluating the Tutte 
polynomial for this graph at x — 1, one encounters an expression of the form 0/0 and must 
use L'Hopital's rule to evaluate it. Thus, differentiating numerator and the denominator 
l/(x — 1) from the prefactor and using eqs. (4.2.8) and (4.2.10), we have 



n T (L y ,d) 

E 

d=l ~ j=l 



L v / , , \ i n T {L y ,d) 

T[(G s (L y xm ] FBC y ,PBC y ),x=lM = Y,(- 1 ^ 1 (y- 1 r- J ^ 1 - E ( A ^,^) r 
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Ly n T (L y ,d) 

E ^a^ST 1 T ^ (4-2.12) 

d=0 j=l 

where G s denotes the type of lattice and the right-hand side is evaluated at x — 1. This 
provides one way of understanding the origin of certain factor of m in the exact results that 
we shall give below for reliability polynomials of strip graphs of length m. 
The number nT{L y , d) of A's with a given coefficient c (d) is [28] 

*W = T&£k(£i ^,< d <L v (4.2,3) 
and zero otherwise. The total number N TtLy ^ x of different terms \r,L y ,j in eq. (4.1.1) is [31] 

Ly 

N T , Ly , x = J2n T (L y ,d) . (4.2.14) 

d=0 

For cyclic and Mobius strips of the square, triangular, and honeycomb lattices, we calculated 
this to be [31, 32] 

N W = . (4.2.15) 

For cyclic and Mobius strips of regular lattices G s with arbitrary L y , eq. (4.2.13) shows 
that there is a unique \t,Gs,l v 4 corresponding to the coefficient of highest degree d, 
namely d = L y . We have found [28, 29, 32] that this term is 

X T ,G S ,Ly,d=Ly = 1 • (4.2.16) 

As indicated, since this eigenvalue is unique, it is not necessary to append a third index, as 
with the other A's, and we avoid this for simplicity. For each of the cyclic strips of regular 
lattices G s considered here, one can derive a general formula for aa s ,L y ,L y from (3.2) in 
conjunction with (4.2.16), viz., 

0CG„L y ,L y = (l-p) PG > (4.2.17) 

where 

\E\ - \V\ 

p Gs = 1 1 1 1 . (4.2.18) 
m 

For example, for the cyclic and Mobius strips of the square (sq) lattice of width L y vertices 
and length L x = m vertices, \V\ = L y m and \E\ — (2L y — l)m, so that p sq — L y — 1; for the 
corresponding strips of the triangular (t) lattice, |V| is the same and \E\ — (3L y — 2)m, so 
that p t = 2(L y — 1). For the cyclic strip of the honeycomb (he) lattice of width L y vertices 
and length L x = 2m vertices, we have \ V\ = 2L y m and \E\ — (?>L y — l)m so that ph c = L y — l. 
We next have 
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Theorem 4.2 Consider the reliability polynomial for cyclic strips of the square, triangu- 
lar, or honeycomb lattices, denoted generically G s , of width L y and arbitrary length L x . The 
total number of different terms in (J^.1.2) satisfies the inequality Nn i a a ,L y , a < (^/^)^T,G s ,L y ,x 
and hence 

N R , Gs , Ly , a <l^^) • (4.2.19) 



Proof The proof makes use of the expression (3.2) of the reliability polynomial R(G m , p) 
as a special case of the Tutte polynomial T(G m , x, y) for x = 1 and y — 1/(1 — p), together 
with our earlier result (4.2.11) for the Tutte polynomial of a cyclic strip graph of the square, 
triangular, or honeycomb lattice. Because the prefactor l/(x — 1) in (4.2.11) is singular at 
x = 1 where one evaluates the Tutte polynomial, it is necessary that 

= (4.2.20) 

x=l 

This vanishing condition holds for arbitrary m, which means that it implies pairwise relations 
among various terms \T,G 3 ,L y ,d,j- The only such relations that can yield an overall factor of 
(x — 1) are relations of the form \T,G s ,L y ,d,j = ^T,G a ,L y ,d-e,j' for 1 < d < L y and odd £ in the 
range 1 < £ < d for a given d, since these yield expressions of the form 

c (d) (AT, Gs ,L H ^) m + c ( ^ ) (AT, Gs ,L H ,,-, J O m = {c {d) + C id - £) ){\T, GsM ) m (4-2.21) 

Using (4.2.8), one sees that since £ is odd, the combination + c^^ vanishes at q = 
and hence has a factor of q. Using eq. (3.10), one sees that this is sufficient to get an 
overall factor of x — 1 so as to yield the vanishing in eq. (4.2.20) at x = 1 or equivalently 
to cancel the prefactor l/(x — 1). Because of the pairwise equalities connecting each of the 
A's when evaluated at x — 1, the total number of these terms in the reliability polynomial is 
reduced to half of the number for the Tutte polynomial, eq. (4.2.14). Further equalities are, 
in principle, possible and would further reduce the number of different A's. This completes 
the proof of the theorem. □ 

Two remarks are in order here. First, we note that in all of our calculations, we find 
£ — 1; that is, the pairwise equalities among the A's occur for adjacent values of d. Hence, 
the linear combination in eq. (4.2.21) is + c^ d ~ e ^ = n^ d \ where was given in (4.2.4). 
Second, in all of our calculations, we have observed that the inequality (4.2.19) is realized 
as an equality. Let us illustrate how the pairwise equalities work in some special cases. For 
L y = 2, we have ut(2, 1) = 3 and 71^(2, 0) = 2. At x — 1, the single \t,2,2 is equal to one of 
the three \t,2,ij, which we denote At,2,i,3- The other two of these three At,2,i,/s are equal to 
the two respective \t,2,o,/s'- ^t,2,ij = ^t,2,oj for j = 1,2. For L y = 3, we have nr(3, 3) = 1, 
nr(3,2) = 5, n T (3, 1) = 9, and tit(3, 0) = 5. At x — 1, At,3,3 is equal to one of the five 
At,3,2j's while the other four are equal to four of the nine At,3,i,/s, and finally, the remaining 
five of these At,3,ij's are equal, respectively, to the five Xt,3,o/ s - 



Uy 

E 

d=0 



n T (L y ,d) 



T,G S ,Ly,dj) 
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4.3 Cyclic Strips of the Square Lattice with Self-Dual Boundary 
Conditions 

We have previously calculated the Tutte polynomial for cyclic strips of the square lattice 
with self-dual boundary conditions [35, 36, 37]. These strip graphs have (i) a fixed transverse 
width L y , (ii) arbitrarily great length L x , (iii) periodic longitudinal boundary conditions, 
and (iv) are such that each vertex on one side of the strip, which we take to be the upper 
side (with the strip oriented so that the longitudinal, x direction is horizontal) is joined by 
edges to a single external vertex. A strip graph of this type will be denoted generically 
as Gfl (where the subscript D refers to the self-duality) and, when its size is indicated, as 
Gd(L v x L x ). For cyclic strips with self-dual boundary conditions we determined the general 
structure of the Tutte polynomial in [36] . We showed that the coefficients were precisely the 
polynomials. The general form of the Tutte polynomial is [36, 37] 

L y +l n T (G D ,Ly,d) 

T(G D [L y x m,cyc.],x,y) = £ «M £ i A/,,,,.,,,,,)"' (4.3.1) 

d=l j=l 

where R^ was defined in eq. (4.2.5) and 

ML„d) = ^^( T 2L 't},) (4.3.2) 



L y + d + 1 \L y -d+1 
so that for the total number of terms [36], 

AW., = (^l) ' <«- 3 > 

Note that, in contrast to the situation with the cyclic strips, because the Tutte polynomial 
does not contain any prefactor of l/(x — 1), there is no need to use L'Hopital's rule in 
evaluating this polynomial to get the reliability polynomial. Related to this, the coefficients 
are R^ d \ i.e., do not contain the factor of q that the CO efficients do. In this evaluation at 
x — 1, i.e., q — 0, the R^ simply reduce to the constants (— l) d+1 d. The reliability polynomial 
for the cyclic self-dual strip, R{Gn[L y x m],p], thus contrasts with the reliability polynomial 
for the cyclic strip graphs in not having any terms proportional to to, a consequence of the 
fact that one did not have to use L'Hopital's rule in evaluating the corresponding Tutte 
polynomial at x — 1. 
Our next result is 

Theorem 4.3 The reliability polynomial of the self-dual strip graph Gr>\L y x L x ] is (with 
L x = m) 

Ly+l n T (G D ,L y ,d) 

R(G D [L y x m],p) = p L * m E ^G D ,L yAl T (4-3.4) 

d=l j=l 

where 

aG D ,L y 4,j = (1 ~P) Ly ^T,G D ,L y Ai\x=l,y=l/{l-p) • (4.3.5) 
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Proof This result follows directly from (4.3.1) and (3.2), using (4.2.7) and the fact 
that the number of vertices and edges for a graph in this family are \V\ = L y L x + 1 and 
\E\ — 1L y L x . D 

From our explicit calculation of R(Go[L y x m],p) for several values of the width L y , we 
find that there is no reduction in the total number of A's when one carries out the evaluation 
at x — 1 for eq. (4.3.4). Hence, for the cases that we have studied, this total number of A's 
in the reliability polynomial for the self-dual cyclic strips of the square lattice is the same 
as the number for the Tutte polynomial (Potts model partition function) for these graphs. 
Specifically, this is N r ^g d ,l v ,\ = 3, 10, 35 for L y = 1, 2, 3, respectively. 

4.4 Strips with Free Longitudinal Boundary Conditions 

As before for chromatic and Tutte polynomials, it is convenient to express the results in terms 
of a generating function. For a strip graph of type G s (where this includes the specification 
of the transverse boundary conditions) and length m longitudinal edges, we have 

oo 

T(G s ,p,z) = J2R((Gs) m , P )z m (4.4.1) 

m=0 

where 

D(G s ,p,z) 

The numerator and denominator are rational functions of z and p of the form 

dcg z (A0 

Af(G s ,p,z)= A G*jz j ( 4 - 4 -3) 

j=i 

dcgJV) 

V(G s ,p,z) = l+ b Gs,jZ J ■ (4-4.4) 

3=1 

The property (4.1.5) is manifested in the relations 

b Gs ,i = -1 , b G , d = , 2 < j < N Gs , a for p = 1 (4.4.5) 
For the strips of the square and triangular lattices with free longitudinal boundary conditions, 

A Gsfi = R(T Ly ,p)=p L y~ 1 for G s = sq,t (4.4.6) 



4.5 Accumulation Set of Zeros of R(G,p) 

Concerning the continuous accumulation set B of the zeros of R(G,p) for cyclic and Mobius 
lattice strip graphs in the limit L x — > oo, we note 

Theorem 4.4 Consider the lattice strip graphs G[L y , L x ,cyc] and G[L y , L x , Mb}. The 
continuous accumulation set B defined in the limit L x — > oo is the same for both of these. 
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Proof This is a corollary of earlier theorems that we have proved for the accumulation 
set of the zeros of the Tutte polynomial or equivalent Potts model partition function, stat- 
ing that these are the same for a given strip with cyclic or Mobius longitudinal boundary 
conditions. □ 

We next discuss the construction of families of graphs G for which the reliability polyno- 
mial R(G,p) has zeros that extend to the circle \p — 1| = 1. A useful lemma can be stated 
as follows [5, 15]: Let G = (V,E) be a connected graph with no multiple edges. We also 
assume that G has no loops; this assumption incurs no loss of generality since loops do not 
affect the reliability polynomial. Now let G t be the graph obtained from G by replacing each 
edge by £ edges connecting the same pair of vertices. Then 

R(G e ,p)=R(G,p') (4.5.1) 

where 

p' = !-{!- p) 1 . (4.5.2) 

In addition to the proofs in [5, 15], this result can be understood easily from a physics 
viewpoint, using the expression of the reliability polynomial as a special case of the Potts 
model partition function. Replacing every edge in G by £ edges connecting the same vertices 
has the effect of replacing the spin-spin coupling J by J' = £J in the Potts model Hamiltonian 
H, and hence K — > K' = £K, so that, since y — e K (cf. eq. (3.9)), 

y^y> = y\ i e ., U -> ?/ = (tf + l)' - 1 (4.5.3) 

and 

Z(G e ,q,v) = Z(G,q,v') (4.5.4) 

The result (4.5.1) then follows via (3.11) (or (3.2)). 

As noted in [5, 15], using this result, one can generate families of graphs, namely those 
with replicated edges, that extend to \p — 1| = 1. The relation (4.5.4) is equivalent to 

V — r J T(G,x',y') (4.5.5) 

where y' was given in (4.5.3) and 

x' = 1 + fr-^- 1 ) (4.5.6) 

Now let G = (V, E) be a simple graph (i.e. a graph without loops or multiple edges) that is 
connected and has at least one edge. As before, the exclusion of loops does not incur any 
loss of generality since loops have no effect on the reliability polynomial. Now let Ge be 
the graph obtained from G by replacing each edge by £ edges connecting the same pair of 
vertices. Then the zeros of R(Ge,p) in the complex p plane include the following: 

p=l-e™ , k = 0, !,...,£-! (4.5.7) 
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and these satisfy \p — 1| = 1. This follows by noting that in general, for the specified G with 
at least one edge, R(G,p) contains a factor pl^l -1 and hence Gg contains a factor of 
Therefore R{Gi,p) has at least one zero at p' = 0. Recalling the relation (4.5.2) and solving 
for the £ roots of this equation in the p plane, we obtain (4.5.7), which satifies \p — 1] = 1. 

Let us give an example. Perhaps the simplest is to consider the graph G = T2. Replacing 
the single edge in this graph by £ edges yields a graph commonly denoted the "thick link" 
or "fat link", TL e . This is the (planar) dual to the circuit graph Q; i.e., TL e = C\. From 
the elementary result R(T 2 ,p) = p and our theorem (4.5.1), one gets R(FL e ,p) = p', which 
has zeros at the points (4.5.7). The location of these zeros for £ = (a)2, (6)3, and(c)4 are (a) 
p = 0, 2; (b) p — and p = 3/2 ± (i/2)y / 3; (c) p — 0,2, and 1 ± i. In the limit £ — > 00, one 
then obtains an accumulation set B that is the circle \p — 1| = 1. Note that in this example, 
although the number of edges goes to infinity, the number of vertices is fixed at 2. One can 
also consider families of graphs with a subset of the vertices replicated. The family of L y = 2 
strips of the square lattice with cylindrical or torus boundary conditions are of this type, 
with double vertical and single horizontal edges. These families will be analyzed below. 

5 Calculations for Specific Families of Graphs 

5.1 L y = 1 Strips of the Square Lattice 

The reliability function for the free L y = 1 strip is a special case of the result for a tree graph 
T rn with m vertices, namely R = p m ~ l . For the case of cyclic boundary conditions, the strip 
is the circuit graph, C m . By an elementary application of the deletion-contraction theorem, 
one has 

R{C rn ,p)=p m - l [m{l-p)+p\ . (5.1.1) 

Hence, 

r = p (5.1.2) 

for the infinite-length limits of both the free and cyclic L y = 1 strips. The polynomial 
R(C m , p) has a zero of multiplicity m — 1 at p = and another zero at p = mj [m — 1), which 
monotonically approaches p = 1 from above as m — > 00. The r function (5.1.2) for these 
L y = 1 strips has the special property that dr/dp = 1. As will be seen, the r functions that 
we shall calculate for other strips have various nonzero slopes at p = but, in contrast to 
(5.1.2), they have the derivative dr/dp — at p — 1. 

5.2 L y = 2 Strip of the Square Lattice with Free Boundary Con- 
ditions 

The generating function for this strip is of the form (4.4.2) with M = p and 

V{p,z) = l-p 2 (4-3p)z + p*(l-p)z 2 

= (1 - p 2 a sq , 2 ,o,iz)(l - p 2 a sq , 2 ,o,2z) (5.2.1) 
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where 



with 



i? 2 = 12 - 20p + 9p 2 . 



J = 1,2 



(5.2.2) 
(5.2.3) 



The subscripts in the notation make reference to the structure of the a's for the cyclic strip; 
see below. We observe that R(sq[2 x m],p) has the factor (4 — 3p) if m > 1 is odd. For the 
cases that we have examined, besides the real root at p — 0, this reliability polynomial has 
a single real root at p — 4/3 if and only if m > 1 is odd. 

The dominant term in the physical interval p G [0, 1] is a sgj 2,o,i) an d hence the asymptotic 
reliability per vertex for L x — > oo is 



r(sg[2 x oo, free\,p) = p(a sq ^, ,i) 1/2 



V 

V2 



4 - 3p + V 12 - 20p + V 



1/2 



This has the derivatives 



and 



dr 
dp 



p=0 



1 + vg 

v/2 



1.93185. 



dr 
dp 



. 



(5.2.4) 



(5.2.5) 



(5.2.6) 



The locus £> for the L y — 2 free strips of the square lattice, shown in Fig. 1, is formed 
as the continuous accumulation set of zeros of the reliability polynomial in the limit of 
infinite strip length, m — > oo and is the solution of the equation expressing the degeneracy 
of magnitudes 1 0^,2,0,1 1 — l a sq,2,o,2|- This is an arc of a circle with radius 1/3 centered at 1, 



1 



p = i+ e * ee[-e sq2 ,e sq2 ] 



where 



9 sq2 = arctan(2 v / 2 ) ~ 70.53° . 
This circular arc crosses the real p axis at a single point, which is thus p c , namely 



Pc 



for sq, FBCy, 2 x 00 . 



(5.2.7) 



(5.2.8) 



(5.2.9) 



This applies here for the case of free longitudinal boundary conditions, and we shall show 
below that it also applies for the case of periodic or twisted periodic longitudinal boundary 
conditions. We find that this property that the value of p c for a strip of a given type of 
lattice and a given set of transverse boundary conditions is independent of the longitudinal 
boundary conditions also holds for other strips, and we indicate this when listing these values 
of p c (or (p c )eff) by displaying only the transverse boundary condition (FBC y in eq. (5.2.9)). 
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Im(p)0 




-0.5 



Figure 1: Singular locus B for the L x — > oo limit of sg(2 x L x ) for either free, periodic, or twisted periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
cyclic strip with L x = 30 (i.e., n = 60). 
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The arc has endpoints at the values of p where the the expression 12 — 20p + 9p 2 in the square 
roots in a sqj2 ,o,j vanishes so that a agj 2,o,i — a sq,2,o,2- These endpoints are 

p = 1 + - e ±ie °« 2 = -(5±V2i)~ 1.11111 ± 0.31427i . (5.2.10) 
3 9 

Recall that two roots of an algebraic equation are equal where the discriminant vanishes; 
for the equation a 2 — (4 — 3p)a + 1 — p = that yields a S q,2,o,j, j = 1,2 as its roots, this 
discriminant is 12 — 20p + 9p 2 . 



5.3 L y = 3 Strip of the Square Lattice with Free Boundary Con- 
ditions 

The generating function for this strip follows from our calculation of Z(G, q, v) for this strip 
[33] and has the form (4.4.2) with 



and 



J\f = p 2 1 + p\l - p)z - p 6 (l - pfz 2 



(5.3.1) 



where 



V(p,z) = 1 + ^6 



sq3f,j A 



3=1 



24p + 10p 2 

2 



&sg3/,l = -P 3 (15 

fe sg3 /,2 = p\l - p)(32 - 66p + 46p z - llp J 
& sg3 /,3 = V(l - P) 3 (15 - 20p + 7p 2 ) 



P"(1-P) E 



(5.3.2) 

(5.3.3) 
(5.3.4) 
(5.3.5) 
(5.3.6) 



In Fig. 2 we show the locus B for the L y = 3 strips of the square lattice with any 
longitudinal boundary conditions, together with zeros for the cyclic strip of this width with 
L x = 20. This locus consists of two complex-conjugate arcs which do not cross the real 
axis. From the arc endpoints closest to the real axis, one can infer the effective quantity 
{p c )eff — 1-335. The locus is again concave to the left and roughly centered about the point 
p=l. 



5.4 L y = 2 Cyclic and Mobius Strips of the Square Lattice 

For the reader's convenience, we shall give some details of our calculation for this family 
of graphs. The Tutte polynomials for these families of graphs were computed in [28] and 
have Nt,2,\ = 6, in accordance with the general formula (4.2.15). There are nr(2,0) = 2 
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Im(p)0 




-0.5 



Figure 2: Singular locus B for the L x — > oo limit of sg(3 x L x ) for either free, periodic, or twisted periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
cyclic strip of this width with L x — 20 (i.e., n — 60). 
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sq,2,2 



= 1 



terms with coefficient c^°\ denoted At,s ? ,2,oj, j = 1)2; nr(2, 1) = 3 terms with coefficient 
S l \ denoted \T, S q,2,i,j, and the unique term with coefficient S 2 \ which term is At, 
as a special case of the general result (4.2.16). Explicitly, 

j 2 n T (2,d) 

T(sq[L y = 2, L x = m, cyc},x, y) = —— c (d) (At,,,^-) 

Ob X 

d=0 i=l 

"2 3 

>>/..,.2.n.,r + C (1) X)( WW)™ + C (2) (AT, S9 ,2, 2 ) m 



where riT(L y ,d) was recalled from [31] in eq. (4.2.13) and 
1 



(5.4.1) 



A 



T,sq,2,0,j 



1 + y + x + x 2 ) ± (y 2 + 2y(l + x - x 2 ) + (x 2 + x + 1) 



1/2 



Ar,sq,2,l,j 



x + y + 2 ± ( (x - yf + 4(x + y + f 



1/2 



with j — 1,2 corresponding to ±, and 



Ar, sg ,2,i,3 - ^ • 



(5.4.2) 
(5.4.3) 

(5.4.4) 



We have discussed above how equalities occur between certain \t,g,l v 4,Y s when evaluated 
at x — 1. We show this explicitly here. First, we have 

1 



^T, S q,2,o,j = Ar, a9 ,2,ij = - [y + 3 ± vV + 2y + 9 



for x = 1 (5.4.5) 
where j — 1,2 refer correspond to ± on the right-hand side, respectively. Furthermore, 

Ar, sg ,2,i,3 = Ar>g,2,2 = 1 for x = 1 . (5.4.6) 

Hence, N R>sq:2> \ = (l/2)N T , sq ,2,\ = 3, in accordance with the inequality (4.2.19) (realized as 
an equality) in our theorem above. Using L'Hopital's rule to evaluate the Tutte polynomial 
(5.4.1) at x = 1, we obtain 



R(sq[L y = 2,m,cyc},p)=p 2m 
+mp -\l- p ) 2 !(a sq , 2Al r- 1 



a sq ,2fl,i) m + (a sg ,2fi,2) m - 2(a S g, 2 ,2) r 
, , 3(1 -pY 



m—l 



3(1 ~P) 



i a sq,2,2) 



R2 

m—l 



(5.4.7) 



where R 2 was defined above in eq. (5.2.3), the a s?; 2,o,j were given in (5.2.2) above, and 



a.sq,2,2 = 1~P 



(5.4.* 



as a special case of (4.2.17) above. 
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For the L y = 2 Mobius (Mb) strip of the square lattice 



T(sq[L y = 2,L X = m,Mb],x,y) 



x - 1 



y~l(Ar,s g ,2,o,j) T 



■3=1 



6 

+C (1) (-(AT, Sg ,2,l,l) m + E( WW)™) - ( W,2) ? 



i=2 



(5.4.9) 



From this we obtain 



R(sq[L y = 2, to, M6],p) = p 



2m 



(« Sg ,2,0,l) m + («, ?> 2,0,2) m " (a Sg ,2, 2 ) r 



(a S g,2,0,2) 



m— 1 



3(1 -p) 



sq,2,2) 



\m—l 



(5.4.10) 



For both of these L y = 2 cyclic and Mobius strips of the square lattice, the dominant 
term in the physical interval p G [0, 1] is <y sq ,2,o,i and hence the asymptotic reliability per 
vertex, r, is the same as that for the corresponding strip with free longitudinal boundary 
conditions, given by (5.2.4): 



r(sq[2 x oo, eye], p) = r(sq[2 x oo, Mb.], p) 

= r(sq[2 x oo, free],p) 



= r{sq[2 x oo, FBC y ],p) . 



(5.4.11) 



For the infinite-length limit of these strip graphs, and for the respective infinite-length limits 
of other strip graphs to be considered below, we find that the r function is independent of 
the longitudinal boundary conditions. This is indicated in the last line of eq. (5.4.11) and 
analogous equations below by listing only the type of transverse boundary condition. The 
dominant term for large positive p is a S g,2,o,2- Because the two dominant a's in the reliability 
polynomials for these cyclic and Mobius strips are the same as the two a's that enter in the 
reliability polynomial for the corresponding strip with free boundary conditions, it follows 
that the locus B is identical for the L y = 2 strips of the square lattice with free, cyclic, and 
Mobius boundary conditions. This locus is shown in Fig. 1. 

We find that the reliability polynomials R(sq[2 x L x ,cyc.],p) and R(sq[2 x L x , Mb],p) 
have only p = as a real root if L x > 1 is odd, while for even L x > 2, they each have 
an additional real root. In Table 1 we list the values of this additional respective real root 
for the reliability polynomials of these two strips. The results in this table suggest that as 
L x — > oo, these two respective additional real zeros for the cyclic and Mobius strips could 
approach the same limit. The value of this limit could be about 1.3. 
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Table 1: Values of the real root (aside from p = 0) of R(sq[2 x L x ,cyc.],p) and R(sq[2 x L x ,Mb],p) for 
even L x . 



L x 


eye. 

■J 


Mb. 


9 


2 000000 


1 626^38 


4 


1.404813 


1.388652 


6 


1.349507 


1.346559 


8 


1.333333 


1.332635 


10 


1.327103 


1.326919 


12 


1.324413 


1.324362 


14 


1.323229 


1.323214 


16 


1.322761 


1.322757 


18 


1.322658 


1.322657 


20 


1.322749 


1.322748 



5.5 L y = 2 Cylindrical Strips of the Square Lattice 

This is the minimal-width strip of the square lattice with cylindrical boundary conditions 
and involves double vertical edges. The results for the reliability polynomial are conveniently 
expressed in terms of the generating function (4.4.1) and (4.4.2). We find 

M = R(C 2 ,p) = p(2 - p) (5.5.1) 



V = l-p 2 (Q-8p + 3p 2 )z + (l-p) 2 p 4 z 2 



= Yli 1 -P 2 a sqc yl2,jz) 

i=i 



where 



a sqcyl2,j 



Q-8p + 3p 2 ± y/(2 - p)(4 - 3p)(4 - Qp + 3p 2 ) 



(5.5.2) 



(5.5.3) 



We observe that R(sq[2 x L x ,cyl.],p) always has the factor (2 — p); this is associated with 
the fact that this strip has double vertical edges. For odd m > 1, i.e., even L x > 2, 
R(sq[2 x L x ,cyl.],p) also has the factor (3p 2 — 8p + 6), with roots (l/3)(4 ± y/2i). For 
m > 4, we observe that R(sq[2 x L x ,cyl.],p) has, besides the real zeros at p = 0,2, a set 
of real zeros in the interval (1,2]. Among other features, we note that as m increases, (i) 
the number of real zeros in this interval increases; (ii) the smallest zero in this interval 
decreases monotonically toward 4/3; (iii) the largest zero increases monotonically toward 
2 as m increases, and (iv) the density of zeros is largest in the regions slightly above 4/3 
and slightly below 2. A plot of these zeros is shown in Fig. 3; as is evident from this plot, 
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as m — > oo, these real zeros away from p = accumulate to form the line segment on B 
extending between p = 4/3 and p = 2. 

In the physical interval p G [0, 1] the dominant term is cx S qcyi2,i, so the asymptotic relia- 
bility per vertex is 

r(sq[2 x oo, PBC y ],p) = p(a sqcyl2)1 ) 1/2 



V 

V2 



, -I 1/2 

6 - 8p + 3p 2 + y/(2 - p)(4 - 3p)(4 - 6p + 3p 2 ) (5.5.4) 



This has derivatives 



dp 



= 1 + ^ = 2.41421... (5.5.5) 

p=0 

and dr /dp — at p = 1. 

The locus £>, shown in Fig. 3, is given by the union of an arc of a circle with a line 
segment on the real axis: 

B: {p=l + ^=e w , 9 e [-tt/2, tt/2]} U < p < 2} . (5.5.6) 

Hence, 

p c = 2 for sq, PBC y , 2 x oo (5.5.7) 

The circular arc intersects the real axis and line segment at p = 1 + ~ 1.577. The 
endpoints of the arc are located at 

Psqce,P* sqce = l±^ (5-5.8) 

These points, together with the endpoints of the line segment, p = 4/3 and p = 2 are the 
zeros of the expression in the square root in a sqcy i2j in eq. (5.5.3), which expression is the 
discriminant of the equation a 2 — (6 — 8p + 3p 2 )a + (1 — p) 2 = whose roots are these a's. 

5.6 L y = 2 Torus and Klein Bottle Strips of the Square Lattice 

Either using a direct calculation or as a special case of our previous calculation of the Tutte 
polynomial [32], using (3.2), we find 



R(sq[L y = 2, to, torus] , p) = p 



2 m 



(&sqtor2,l) m + (® S qtor2,2) m ~ sqtort^Y 



- mp -. (1 _ pf (°^">"-' 2p _ 3 + ( P -2)(6p»-13 P + 8) 

} 2 ^ r V(P- 2 )( 3 P- 4 )( 3 P 2 -6p + 4) 

( 2p _ 3 _ <^;f^-fr^ 4 ) + (i - ,)( w~ 

a/(p - 2) (3p - 4) (3p 2 - 6p + 4/ 



(5.6.1) 
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Figure 3: Singular locus B for the L x — > oo limit of sqcyl(2 x L x ) for either free or (twisted) periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
torus strip with L x = 30 (i.e., n — 60). 
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where 

®sqtor2,j = &sqcyl2,j ] ' = 1, 2 (5.6.2) 

with a sqcy i2j, j = 1,2, given above in (5.5.3) and, 

a sq tor2,3 = (1 -P? • (5.6.3) 
For the L y = 2 Klein bottle (Kb.) strip of the square lattice 



R(sq[L y = 2,m,Kb.],p) = p 



2m 



{&sqtor2,l) m + {°^sqtor2,2) m ~ sqtor2 fl)™ 



_ mp -. (1 _ p f (^J- 2p _ 3 + ( P -2)(6^-13 P + 8) 

j 2 ^ v/(P-2)(3p-4)(3p 2 -6p+4) 

V(p-2)(3p-4)(3p'-6p + 4)i *"*»• 



m— 1 



(5.6.4) 



In the physical region p G [0, 1], the dominant term is a sq t or 2,i for both of these families of 
strips, and hence r(sq[2 x oo, torus/Kb.}, p) = r(sq[2 x oo, cyl.],p). This equality has already 
been incorporated in our notation r(sq[2 x oo, PBC y ],p) in eq. (5.5.4). This is another 
example of the feature that for a given type of strip graph G s with a given width and set of 
transverse boundary conditions, r is independent of the longitudinal boundary conditions. 

Similarly, since the locus B is determined by the equality in magnitude of the two dom- 
inant eigenvalues a sq t or 2,i and a sqt0 r2,2, and since these are the same for the corresponding 
strip with cylindrical boundary conditions, it follows that this locus is the same as the locus 
for the L y = 2 cylindrical strip of the square lattice, shown in Fig. 3. 

5.7 L y = 3 Cylindrical Strips of the Square Lattice 

Again we express the results in terms of the generating function for the reliability polynomial, 
(4.4.1) and (4.4.2). We find 

J\f = R(C 3 ,p) + p 5 (3 - p)(l - p) 2 z = p 2 [3 -2p + p 3 (3 - p)(l -p) 2 z] (5.7.1) 

where RiC^.p) = p 2 {3 — 2p), and 

V = 1 - p 3 (24 - 56p + 46p 2 - \3p 3 )z + p 6 (l - p) 2 (24 - 4Qp + 30p 2 - 7p 3 )z 2 - p 9 (l - pfz 3 

3 

The largest of the roots of the associated equation 

£ 3 - (24 - 56p + A6p 2 - I3p 3 )£ 2 
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Table 2: Values of the real root (aside from p = 0) of R(sq[3 x oo, cyl.],p) for odd L, 



T 


root 


1 


1.500000 


3 


-i A ,1 n n n f—r 

1.440337 


5 


1.422909 


7 
( 


1/11 fiOQ7 


9 


1.412887 


11 


1.410809 


13 


1.409409 


15 


1.408403 


17 


1.407644 


19 


1.407051 


21 


1.406575 



+ (1 - p) 2 (24 - 46p + 30p 2 - 7p 3 )£ - (1 - pf = (5.7.3) 

is the dominant a sqcy i 3)max , so that 

r(sg[3 x co,PSC y ],p) = p(a sgcy i 3 , maa; ) 1/3 . (5.7.4) 

This function has the property that dr/dp = 2.84207... at p = and dr/dp — at p — 1. 
We find that the reliability polynomial i?(sg[3 x oo, cyl.],p) has only p = as a real root if 
m > 1 is odd, i.e. L x > 2 is even, while for even m > 0, i.e., odd L x > 1, it has an additional 
real root which decreases as m increases and appears to approach a limit of roughly 1.4. In 
Table 2 we list the values of this additional real root. 

In Fig. 4 we show a plot of B for the L x — > oo limit of the L y = 3 strip of the square 
lattice with cylindrical boundary conditions. 

5.8 L y = 4 Cylindrical Strip of the Square Lattice 

As before, we express the results for the L y = 4 strip of the square lattice with cylindrical 
boundary conditions in terms of the numerator and denominator of the generating function. 
In particular, for the denominator, we find 

6 
3 = 1 

6 

= Yli 1 ~ P 4 ®sqcyl4,jZ) (5.8.1) 
3=1 
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Figure 4: Singular locus B for the L x — > oo limit of sqcyl(3 x L x ) for either free or (twisted) periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
toroidal strip with L x = 20 (i.e., n = 60). 
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where 

b S qayi4,i = -p 4 (90 - 293p + 369p 2 - 211p 3 + AQp 4 ) (5.8.2) 
b sqcy i4,2 = p 8 (l - p) 2 (735 - 2977p + 5094p 2 - 4710p 3 + 2481p 4 - 706p 5 + 85p 6 ) (5.8.3) 

b sqcy i4,3 = -P 12 (l - #(1548 - 7518p + 15948p 2 - 19170p 3 

+ 14104p 4 - 6351p 5 + 1621p 6 - 181p 7 ) (5.8.4) 

bsqcyUA = P m (l ~ #(735 - 3177p + 5870p 2 - 5934p 3 + 3456p 4 - 1097p 5 + 148p 6 ) (5.8.5) 
b sqC yu,5 = ~P 20 (1 - P) H (90 - 223p + 208p 2 - 88p 3 + 15p 4 ) (5.8.6) 
b sqcy iA fi = p 24 (l-p) 15 (5.8.7) 

The degree of T> is too high to obtain closed-form algebraic expressions for the ot sq cyi4,jS) 
but they may be computed numerically. Denoting the dominant Oi S qcyi4,j i n the physical 
interval p G [0, 1] as a sqcy i^ max , one has, for the asymptotic reliability per vertex, r(sq[4 x 
oo, PBCy],p) = p^sqcyi^max) 1 ^ ■ In Fig- 5 we show the locus B for the infinite-length limit 
of this strip. This locus includes two complex-conjugates arcs, a self-conjugate arc crossing 
the real axis at p ~ 1.364, and a real line segment extending from p = 4/3 to p ~ 1.384, 
which defines the value of p c for this strip. 



5.9 L y = 2 Strip of the Triangular Lattice with Free Boundary 
Conditions 

In addition to studying how reliability polynomials and their asymptotic r functions depend 
on the width and boundary conditions, it is also of interest to explore how they depend on 
the lattice type. For this purpose we have carried out calculations of reliability polynomials 
for strips of the triangular and honeycomb lattice. We begin with the L y = 2 strip of 
the triangular lattice with free boundary conditions, for which we calculate the generating 
function with 

Af(p,z) =p l+p(l-p)z (5.9.1) 

V(p,z) = 1 -p 2 {7- 10p + 4p 2 )^ + j9 4 (l -p)V 

= (l-p 2 a t ^ 1 z)(l-p 2 a t , 2 ,o,2z) (5.9.2) 

a t ,2,oj = \[7- 10p + Ap 2 ± (3 - 2p) ^5 - 8p + 4p 2 ] 

-3-2p + ^-8p + 4p*- for . = 1>2 _ (5Q3) 



where 



28 



Im(p)0 




-0.5 



Figure 5: Singular locus B for the L x — > oo limit of sqcyl(A x L x ) for either free or (twisted) periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
cylindrical strip with L x = 21 (i.e., n — 84). 



29 



We note that R(p, z) for the strip of the triangular lattice of length m has the factor (8 — 
lip + 4p 2 ) if m — 1 mod 3, where our convention is that m = 1 is a single square with 
diagonal (i.e. two triangles), and so forth. 

The dominant term in the physical interval p G [0, 1] is 0^2,0,1, and hence the reliability 
per vertex for L x — > 00 is 



r(t[2x 00, FBC y ],p) = pK 2 ,o,i) 1/2 



This has derivatives 



dr 
dp 



P 

V2 



p=0 



2.61803.. 



and dr /dp — at p — 1. 

The locus B, shown in Fig. 6, is given by an arc of a circle: 

B: ,-l + ie", » 6 [-=,=] 
This is a leftwardly concave circular arc that crosses the real axis at 



1/2 



7 - 10p + 4p 2 ± (3 - 2^5 - 8p + 4p 2 " (5.9.4) 
3 + ^ 



(5.9.5) 



(5.9.6) 



p c = - tri, FBCy, 2 x 00, 



(5.9.7) 



and has endpoints at p — 1 ± i/2, where the factor 5 — 8p + 4p 2 in the discriminant (5 — 
8p + 4p 2 )(3 — 2p) 2 of the equation a 2 — (7 — lOp + 4p 2 )a + (1 — p) 2 = yielding the a's 
vanishes. We shall see that a rather different strip, the self-dual strip of the square lattice 
of width L y — 1, yields the same locus B, although the reliability polynomials are different. 
Interestingly, the effective degree is the same, namely 4, for these two strips, although the 
L y = 2 cyclic strip of the triangular lattice is a A-regular graph, with uniform vertex degree 
A = 4, whereas the self-dual strip of the square lattice with L y — 1 has two quite different 
kinds of vertex degrees - 3 for the outer rim and m for the central vertex connected to each 
vertex on the rim by edges. However, it is not, in general, true that if two families of graphs 
have the same effective vertex degree then the resultant loci B are the same; for example, 
each of the self-dual cyclic strips of the square lattice of width L y has the same effective 
vertex degree in the limit L x — > 00, but they have different loci B. 



5.10 L y = 2 Cyclic Strip of the Triangular Lattice 

We next consider the L y = 2 cyclic strip of the triangular lattice. Either using a direct 
calculation or as a special case of our previous calculation of the Tutte polynomial [29], 
using (3.2), we find 



R(t[L y = 2,m,cyc},p) = p 



2m 



K 2 ,o,i) m +(« W ) m -2(« w )' 



30 



1 



0.5- 



lm(p)0"- 



-0.5- 




0.5 1 1.5 2 

Re(p) 



Figure 6: Singular locus B for the L x — > oo limit of tri(2 x L x ) for cither free or (twisted) periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
cyclic strip with L x = 30 (i.e., n = 60). 
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+mp-\5 - 3 P )-\1 - pf I K^r^l - p) (l - Ap + ^=== 

+ K 2A2 r- i (i-p) ^-^- { ^^==) -2(p 2 -3 P+2 )K 2i2 r- 1 

(5.10.1) 

where «i,2,o,i were given above for the free strip and 

04,2,2 = (1 - P? ■ (5.10.2) 

Because the dominant a's consist of at,2,o,j, j — 1) 2, which are the same for the free, cyclic, 
and Mobius strips, it follows that the locus B is the same for the L x — > oo limits of these 
three strips. 

5.11 L y = 2 Free Strips of the Honeycomb Lattice 

For the L y = 2 strip of the honeycomb lattice with free boundary conditions we calculate a 
generating function with 

M = p (5.11.1) 



V = l-p A (Q-5p)z+(l-p)p s z 2 
= (1 - p 4 a hc , 2 ,o,iz)(l - p 4 a hc , 2 ,o,2z) 



where 



with 



1 r 



6 - hp ± V# 2 



(5.11.2) 
(5.11.3) 



H 2 = 32 - 56p + 25p 2 . (5.11.4) 

In the physical interval p € [0, 1] the dominant term is a^ C) 2,o,i, so the asymptotic reliability 
per vertex is 



r(hc[2xoo,FBC y ],p) = p(a hc>2 , ,i) 1/4 



P 



Q-5p+ ^32 - 5Qp + 25p 2 



This has derivatives 



and dr /dp — at p — 1. 



dp 



2V4 

1 + v 7 ^) 172 = 1.55377... 



1/4 



(5.11.5) 
(5.11.6) 



p=0 



32 



Im(p)0 




-0.5 



Figure 7: Singular locus B for the L x — > oo limit of /ic(2 x L x ) for cither free or (twisted) periodic 
longitudinal boundary conditions. For comparison, zeros of the reliability polynomial are shown for the 
cyclic strip with L x = 30 (i.e., n = 120). 
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The locus £>, shown in Fig. 7, is given by an arc of a circle: 



B: p = l + -e w , 9e[-9 hce ,9 hce } 
o 



where 



9 hce = arctan(4/3) ~ 53.13° 
This is a leftwardly concave arc that crosses the real axis at 

6 

p c = - for he, FBCy, 2 x oo, 



and has endpoints at 
where the discriminant H2 vanishes. 



Phce,P hc e = 7^( 7±i ) 



4 

25' 



(5.11.7) 
(5.11.8) 

(5.11.9) 
(5.11.10) 



5.12 L y = 2 Cyclic and Mobius Strips of the Honeycomb Lattice 

For the L y = 2 cyclic strip of the honeycomb lattice, either using a direct calculation or as 
a special case of our previous calculation of the Tutte polynomial [32], using (3.2), we find 



R(hc[L y = 2,m,cyc],p) = p 4m 
+m p- 1 (l-pfS(a hc , 2fl , 1 ) m - 1 (3 + 



(«?ic,2,0,l) m + ((Xhc,2,0,2) m ~ 2(0^,2,2)' 

16 - 15p N 

7m 



+(a h c2,o,2) m - 1 ( 3 - ) - 2K Cj2 , 2 r- 1 



'Ho 



(5.12.1) 

where ahc,2,o,j, j = 1,2, were given above in (5.11.3) and, as a special case of (4.2.17), we 
have 

0^,2,2 = 1-?. (5.12.2) 
For the L y = 2 Mobius strip of the honeycomb lattice, using the same methods, we obtain 



R(hc[L y = 2,m,Mb],p) =p 



Am 



(«/ lc ,2,0,l) m + (tthc,2,0,2) m — (Ctyic,2,2) 



- p) 2 \ (afcew)™" 1 ( 3 + 



+ (afe c ,2,0,2) 



/ 3 16 ~ 15 P 



2(0^,2,2) 



m— 1 



(5.12.3) 
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Since the locus B is determined by the equality in magnitude of the two dominant eigen- 
values ck/j C) 2,o,i and 0^,2,0,2, and since these are the same for the strip with free and periodic 
or twisted periodic boundary conditions, it follows that this locus is the same as (5.11.7) and 
(5.11.8) for the L y = 2 cyclic or Mobius strip of the honeycomb lattice. We show this locus 
in Fig. 7. 



5.13 L y = 1 Cyclic Strip of the Square Lattice with Self-Dual 
Boundary Conditions 

We have presented some general structural results above for reliability polynomials of cyclic 
strip graphs of the square lattice with self-dual boundary conditions. Recall that these are 
denoted Go, of size L y x L x ; in subscripts we will use the symbol sqdbc. In this and the 
next section we give explicit calculations of the reliability polynomials. Note that Gd[1x L x ] 
is the wheel graph with a central spoke vertex connected to L x outer vertices on the rim. 
Either using a direct calculation or as a special case of our previous calculation of the Tutte 
polynomial [35, 36], using (3.2), we find (with L x = m) 



R(G D [1 x m\,p) = p m (a sqdbcX1<1 ) m + (a sqdbC)1X2 ) m - 2(a sqdbcX2 , 



where 



a sqdbc,l,l,j 



3 - 2p ± ^5 - 8p + Ap 2 ] j = 1,2 



and, consistent with (4.2.17), 



Oi sq dbc,l,2,l — 1 — P ■ 



(5.13.1) 



(5.13.2) 



(5.13.3) 



In accordance with our general results above, the total number of terms is the same as that 
for the Tutte polynomial, N RjGDjLy=1:X = N TtGDtLy=1 ^ = 3. 

In the physical interval p G [0, 1] the dominant term is cx S qdbc,i,i,i so ^ na ^ ^ ne asymptotic 
reliability per vertex is 



r(G D [l x oo],p) = pa sqdbcX1 ,i 

3-2p+ y/l - 8p + 4p 2 



V 



(5.13.4) 



Interestingly, although the graphs are different and the reliability polynomials are different 
for the Gd(1 x m) and tri(2 x m) strips (with either both free or both periodic longitudinal 
boundary conditions), the infinite-length limits of these four families of graphs yield the same 
r function, i.e., (5.13.4) is identical to r(tri[2 x 00, FBC y ],p) given above in (5.9.4). 

The locus B is shown in Fig. 8 and forms arc of a circle centered at p — 1 of radius 1/2 
with endpoints at p — 1 ± (l/2)z, i.e., 



B : p 



0e[- 



7T 7T 

2' 2 J 



(5.13.5) 
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Figure 8: Singular locus B for the L x — > oo limit of G_d(1 x L x ). For comparison, zeros of the reliability 
polynomial are shown for L x = 40. 
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so that 

p c = - for G D (1 x oo) . (5.13.6) 

The endpoints occur at the branch point singularities of the polynomial 5 — 8p + 4p 2 that 
appears in the square root in the terms (5.13.2). The locus (5.13.5) is identical to (5.9.6), 
although the ce's involved are different. 



5.14 L y = 2 Cyclic Strip of the Square Lattice with Self-Dual 
Boundary Conditions 

ForG D [2xm], our general formulas give N RjGdj2 , \ = 10 with n R (G D , 2, 1) = 5, n R (G D) 2, 2) = 
4, and n R (G D ,2,Z) = I. We find that R(G D [2 x m],p) has the form (4.3.4), 

5 4 

R(G D [2 x m],p) = p 2m [X>*2,ij) m - 2 ^(a sg(ibc , 2 , 2 ,,) m + 3(a S(?d6c , 2 , 3 ,i) m ] (5.14.1) 

j"=i j"=i 

where 

<x aq dbc,2,3,i = (1 - P? (5.14.2) 
and the «sgdb c ,2,2,j are the roots of the degree-4 equation 

£ 4 + a dbl £ 3 + a db2 £ 2 + a db3 £ + a dM = (5.14.3) 

with 

a dbl = -(l-p)(7-5p) (5.14.4) 

a db2 = (l-p) 2 (l3-17p + 5p 2 ) (5.14.5) 

a db3 = -(l-p) 4 (7-5p) (5.14.6) 

a dM = (l-pf ■ (5.14.7) 
In eq. (5.14.1) the ot S qdbc,2,i,j are the roots of the degree-5 equation 

e + b dbl e + b db2 e + b db3 e + « + b db5 = o (5.14.8) 

where 

b dbl = -(12- 18p + 7p 2 ) (5.14.9) 

b db2 = (1 - p)(36 - 67p + Alp 2 - 8p 3 ) (5.14.10) 

& d 63 = -(1 -p) 2 (36 - 75p + 53p 2 - 13p 3 ) (5.14.11) 

b dM = 2(1 -p) 4 (3 - 2p)(2 -p) (5.14.12) 

^ 5 = -(l-p) 6 ■ (5.14.13) 

The locus B is shown in Fig. 9 and consists of arcs, again concave to the left, that almost 
cross, but actually have endpoints near to, the real axis at q ~ 1.4 and have endpoints at 
q ~ 0.7765 ± 0.4302i and q ~ 1.406 ± 0.036795i. 
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Figure 9: Singular locus B for the L x — > oo limit of G_d(2 x L x ). For comparison, zeros of the reliability 
polynomial are shown for L x — 30 (i.e., n — 61). 
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5.15 L y = 3 Cyclic Strip of the Square Lattice with Self-Dual 
Boundary Conditions 

For Gjj[3xm], our general formulas give N R; g d ,3,x = 35 with ur^Gd, 3, 4) = 1, n R (Gc, 3, 3) = 
6, tir(Gd, 3, 2) = 14, and ti r (Gd, 3, 1) = 14. We have calculated the reliability polynomial 
from our earlier calculation of the Tutte polynomial for this strip [36, 37]. Since the results 
for the a's are rather complicated, we do not list them in detail here (some examples are 
given in the appendix) but instead concentrate on the locus B. This is shown in Fig. 10. 



Im(p)0 




Figure 10: Singular locus B for the L x — ► oo limit of Gd(3 x L x ). For comparison, zeros of the reliability 
polynomial are shown for L x — 20 (i.e., n — 61). 



5.16 Some Families of Graphs with Multiple Edges 

Consider a connected graph G = (V,E). We take this graph to be loopless; this incurs no 
loss of generality since loops do not affect the reliability polynomial. Now consider a graph 
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H obtained from G by adding another edge. Clearly, for the physical range p e [0, 1], the 
inequality R(H,p) > R(G,p) holds, with equality only at p — 1. That is, adding further 
communications link(s) between two (already connected) nodes in a network increases the 
reliability of the network. A particularly simple modification in the graph is to replace each 
edge joining two vertices by £ edges joining these edges, thereby obtaining what we have 
denoted Gi above. Given any such graph G, whether or not it is a member of a recursive 
family, our theorem above with (4.5.1) and (4.5.2) enables one to calculate R{G^p). Here we 
shall present some illustrations of this. For definiteness, we consider the L y = 2 cyclic strip 
of the square lattice. Using our calculation of the reliability polynomial for the L y = 2 strip 
of the square lattice together with our theorem, we have calculated the reliability polynomial 
for the same strip with all edges replaced by Mold multiple edges joining the same vertices. 

We have calculated the loci B in the L x — > oo limit for these strips. The single point at 
p = 4/3 where the locus B crosses the real p axis for the sq[2 x oo] case, shown in Fig. 1, 
corresponds to £ separate points for the Mold edge replication of sq[2 x oo]. These can be 
calculated using (4.5.1) with (4.5.2), and we obtain 



p = 1 — 



i7r(l + 2fc) 

e I 



k = 0, 1, ...,£- 1 



(5.16.1) 



For example, for i = 2, this yields the two values p — l±i/ \/3, and so forth for higher values 
of £. The arc endpoints can be calculated in the same way. 

It is also of interest to study strip graphs in which a certain subset of the edges are 
replaced by Mold replicated edges joining the same vertices. A simple example of this is 
already provided by the L y = 2 strips of the square lattice with cylindrical, torus, or Klein 
bottle boundary conditions, which have double transverse edges. We have calculated the 
reliability polynomials for the cyclic L y = 2 strips of the square lattice with Mold multiple 
transverse edges. We denote these strip graphs as sq[2 x m,cyc.,v£}. We find 



R(sq[L y = 2,m,cyc,v£\,p) = p 



2m 



+mp-\l-p)^ {{a sqcyM ) m ~ l (VLl(p,£) + VL2(p,£)) 
+ (« sgcyc2l<2 ) m - 1 (VLl(p,£) - VL2(p,£)) - {a sqcycW ) m - 1 } 



(5.16.2) 



where 



(%sqcyc2,v£,j (2p) 



-i 



2-p-(2-3p)(l-pY 



±[(1 - (1 - pY)((2 - pf - (2 - 3p) 2 (l - pY)] 1 ' 2 



, J = 1,2 



(5.16.3) 



&sqcyc2,ve,3 ~ (1 ~ pY 



(5.16.4) 
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and 



VL2(p,£) 



VLl(p,£) 
(2p-l)(3p-2)(l 



(i- p y~\2 P -i) + i 
2p(i- P y- 1 

■p) e - 4 + lOp 



7p 2 + 



2-P 



(1-P) 



2p[(l - (1 - pY){(2 - pf - (2 - 3p) 2 (l - p)^)]V2 



(5.16.5) 
(5.16.6) 



lm(p)0 




Figure 1 1 : Singular locus B for the L x — > oo limit of the L y = 2 strip of the square lattice with I = 3-fold 
multiple transverse edges. For comparison, zeros of the reliability polynomial are shown for this strip with 
cyclic boundary conditions and L x = 20. 



In Figs. 11 and 12 we show the loci B in the limit L x — > oo for the cases I = 3,4. In 
general, B contains arcs protruding inward from points on the circle \p — 1 1 = 1 given by 
p = 1 — exp(2i7ik/£) (cf. eq. (4.5.7)) for < k < i — 1. If and only if £ is even, then one 
of these points, namely the one corresponding to k = £/2, occurs at p = 2, and the part 
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Im(p)0 




Figure 12: Singular locus B for the L x — > oo limit of the Z/ a = 2 strip of the square lattice with £ = 4- fold 
multiple transverse edges. For comparison, zeros of the reliability polynomial are shown for this strip with 
cyclic boundary conditions and L x = 20. 
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of B that connects with this point is a line segment on the real axis extending from p = 2 
downward to an ^-dependent lower limit. 



5.17 L y = 2 Free Strip of the sqd Lattice with Multiple Transverse 
Edges 

In this section we exhibit the first of several recursive families of lattice strip graphs that 
we have found with the property that some zeros of their reliability polynomials and, in 
the limit of infinite length, some portions of their accumulation sets £>, lie outside the disk 
\p — 1|. This family is a one-parameter family generalizing the lowest-order graph found 
by Royle and Sokal. That graph is constructed as follows [19]: start with the complete 
graph K 4 and choose two nonintersecting edges; replace each of these edges by six edges in 
parallel, i.e., joining the same pair of vertices. (Here the complete graph K n is defined as 
the graph with n vertices such that every vertex is connected by one edge to every other 
vertex.) We had previously calculated Tutte and chromatic polynomials for strips composed 
of K4 subgraphs [51, 50, 34]. One of the equivalent ways of defining these strip graphs is to 
start with a strip of the square lattice and add edges connecting diagonally opposite vertices 
of each square, thereby replacing each square by a K4, they were thus denoted strips of the 
sqd lattice, sg^(2 x L x , BC X ), where BC X indicated the longitudinal boundary conditions. In 
a physical context, the Potts model on the sqd lattice is equivalent to the Potts model on 
the square lattice with next-nearest-neighbor spin-spin couplings. The longitudinal (x) and 
transverse (y) directions on this strip are taken to be horizontal and vertical, respectively. 
In particular, we consider a sqd strip with width L y = 2 and arbitrary length L x . We then 
replace each vertical edge by six edges joining the same pair of vertices (leaving the horizontal 
and diagonal edges unchanged). An elementary proof [50] shows that for free longitudinal 
boundary conditions this is a planar graph. We calculate the reliability polynomial and find 
for the numerator and denominator of the generating function 

Af(p,z)=p[A + A 1 zj (5.17.1) 

where 

A = -{p - 2){p 2 - 3p + 3)(p 2 - p + 1) (5.17.2) 

A x = -p 2 (p- l) 6 (3p-4) (5.17.3) 

V(p, z) = 1 + b x z + b 2 z 2 (5.17.4) 

where 

b x = -p 2 (6p 8 - 52p 7 + 200/ - U8p 5 + 6Up 4 - 616p 3 + 391p 2 - 156p + 32) (5.17.5) 

b 2 = -2p\p-2){p-lf (5.17.6) 

Using this exact calculation, one easily verifies that members of this family have reliability 
polynomials with zeros outside the disk \p — 1| < 1. The locus B for the infinite-length limit 
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of this family is shown in Fig. 13. As is evident, B consists of arcs, and the tips of two 
complex-conjugate arcs on the upper and lower left end at the points 

p = 0.327752 ± 0.747464i (5.17.7) 

which have 

[p — 1| = 1.005296 (5.17.8) 

These arc endpoints and their locally neighboring sections of their arcs thus lie outside of 
the disk \p— 1| < 1. Since B is the continuous accumulation set of the zeros of the reliability 
polynomial, it follows that as L x — > oo, infinitely many zeros of the reliability polynomial lie 
outside this disk. Interestingly, although the locus B does extend outside the disk \p— 1| < 1, 
the amount by which it does so is quite small, as is evident from the positions of the arc 
endpoints given in eq. (5.17.8). This property is also true of the other recursive families that 
we have found with zeros of R(G, p) and limiting loci B lying outside the disk \p — 1| < 1. 




Figure 13: Singular locus B for the L x — > oo limit of the L y — 2 strip of the sqd lattice with 6-fold 
multiple transverse edges. For comparison, zeros of the reliability polynomial are shown for this strip with 
free boundary conditions and L x = 21. Horizontal and vertical axes are Re(p) and Im(p) with tick spacing 
0.25, and the circle \p — 1| = 1 is shown with light dots. 

We have also calculated the reliability polynomial for the analogous strip graphs with 
free longitudinal boundary conditions and with each longitudinal, rather than transverse, 
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edge replaced by 6 edges joining the same pair of vertices. We find similar results. A plot 
of the zeros for L x = 9 is shown in Fig. 14. The arcs on B that extend outside the disk 
\p — 1| < 1 are again situated in the upper and lower left, and the endpoints of these arcs 
that extend outside this disk are 

p = 0.4346475 ± 0.8266808? (5.17.9) 

with 

\p- 1| = 1.001511 (5.17.10) 
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Figure 14: Zeros of the reliability polynomial for the L y = 2, L x = 9 strip of the sqd lattice with 6-fold 
multiple longitudinal edges and free boundary conditions. Horizontal and vertical axes are Re(p) and Im(p) 
with tick spacing 0.25, and the circle \p — 1| = 1 is shown with light dots. 



5.18 L y = 2 Free Strip of K± Subgraphs with Multiple Edges 

We consider the strip shown in Fig. 8 of the Appendix of [50] , but with free, instead of cyclic, 
longitudinal boundary conditions. In the notation of that figure, we replace each of the edges 
(labelled as ey = (vi,Vj)) (1,2), (2,3), etc. on the top; (4,5), (5,6), etc. on the bottom, and 
(1,7), (5,8), (2,9), etc. in the middle with a 6-fold replicated edge. This is a family of planar 
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graphs (as is the family in [50] with cyclic boundary conditions). We calculate the reliability 
polynomial and find for the numerator and denominator of the generating function 

N{p,z) =pU + A lZ \ (5.18.1) 

where 

A) = 1 (5.18.2) 
A 1 = /(p-l)(p-2)(/-p+l)(/-3p+3)(2/-14/+42/-70/+70p 2 -42p+13) (5.18.3) 

V(p, z) = 1 + hz + b 2 z 2 (5.18.4) 

where 

&! = p 2 (6p 13 - 86p 12 + 574/ 1 - 2366/° + 6734/ - 14014/ + 22018/ 
-26561/ + 24731/ - 17661/ + 9471/ - 3647/ + 918p - 118) (5.18.5) 
b 2 = /(l - / 3 (4/ 2 - 100/ 1 + 1201/° - 9224/ 9 + 50876/ 8 - 214544/ 7 
+719020/ 6 - 1965184/ 5 + 4459758/ 4 - 8511200/ 3 + 13781747/ 2 - 19045712/ 1 
+22536020/° - 22849416/ + 19815434/ - 14628596/ + 9116492/ 

-4732392/ + 2004288/ - 670184/ + 167360/ - 28048p + 2401) (5.18.6) 

Zeros of the reliability polynomial for the L x = 9 member of this family are shown in Fig. 15. 
The tips of two complex-conjugate arcs, again on the upper and lower left, extend outside 
the disk \p — 1| < 1, ending at the points 

p = 0.4254334 ± 0.8216255i (5.18.7) 

which have 

\p- 1| = 1.002594 (5.18.8) 

Thus, as L x — ■> oo, an infinite number of zeros of the reliability polynomial accumulate to 
form part of the locus B which extends outside the disk \p — 1| < 1. As before, the amount 
by which some zeros and B lie outside the disk is small. 



6 Complete Graphs 

The complete graph K n is the graph with n vertices such that every vertex is connected by 
edges to every other vertex. The family of complete graphs is not a recursive family, and 
provides a contrast to the recursive families on which we have concentrated in this paper. 
The reliability polynomial for K n is [5] 

R{K n ,p) = l-^( n - l \l-py^R{K v p) . (6.1) 
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Figure 15: Zeros of the reliability polynomial for the L y = 2, L x = 9 free strip of K4 subgraphs with 6-fold 
multiple longitudinal edges. Horizontal and vertical axes are Re(p) and Im(p) with tick spacing 0.25, and 
the circle \p— 1| = 1 is shown with light dots. 
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Since this family is not recursive, the reliability polynomial does not have the form (4.1.2) and 
our usual methods for determining an asymptotic accumulation set of zeros as the solution of 
the degeneracy in magnitudes of dominant a terms do not apply. However, we have studied 
the zeros of R(K n ,p) for a wide range of values of n, and we find that these zeros typically 
form an oval-shaped pattern that is concave to the left and is roughly centered around p — 1 , 
somewhat similar to the exact results that we have established for our recursive families of 
graphs without multiple edges. As illustrations, we show in Fig. 16 the patterns of zeros 
for n — 5, 10, and 15. We observe that as n increases, the oval moves outward from p — 1. 
Note that since K n is a A-regular graph with A = n — 1, it follows that the vertex degree 
A — > oo as n — > oo. This feature that A — > oo as n — > oo is thus shared in common by 
the two families K n and the thick links TL n . For the latter family, it is elementary that as 
n — > oo, B is the circle \p — 1| = I. However, our results show that it is not necessary for 
any vertex degree to go to infinity in order for a part of B to satisfy \p — 1| > 1. Our results 
for the family of wheel graphs, Gjj(1 x m) also show that the property that a vertex degree 
goes to infinity as \V\ — > oo is not sufficient for B to have a part with \p — 1| > 1. 

7 Discussion of Results 

In this section we discuss some general features of our results. 

• For a given type of lattice and a given width and choice of transverse boundary condi- 
tions, we find that in the infinite-length limit, the reliability per vertex r is the same 
for any choice of longitudinal boundary conditions. 

• On general grounds, one expects that the higher the connectivity, and, equivalently for 
our present purposes, the higher the effective vertex degree d e ff, the larger the value 
of r for a fixed p. In Fig. 17 we plot r for some strips of the square, triangular, and 
honeycomb lattices. We see that our expectation is confirmed by our exact results. 

• Closely related to this, our results exhibit the property that the derivative dr/dp at 
p = is an increasing function of d e ff. Except for the L y — 1 strips of the square 
lattice, i.e., the line and circuit graphs, we find that dr/dp = at p = 1 for all strips 
that we have studied. In Table 3 we list some of the properties that we have found 
from our calculations. 

• We can also comment on features of the zeros of the reliability polynomials and their 
continuous accumulation sets B. We find that the zeros of the reliability polynomials 
for many strips studied do exhibit the property that \p — 1| < 1, as do their loci B in 
the limit of infinite length. However, motivated by [19], we have calculated reliability 
polynomials for strips such that the zeros of R{G,p) violate the bound \p— 1| < 1 and 
such that, in the limit of infinite length, the continuous accumulation sets of zeros B 
extend outside of the disk \p — 1| < 1. An intriguing result is that the amounts by 
which some individual zeros, and the outer part of the loci B, lie outside this disk are 
small. 
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We find that in some cases, B crosses the real axis and hence defines a p c , but in 
other cases, it does not. However, typically, even for families for which B does not 
cross the real axis, there are arc endpoints on B that lie close to this axis, thereby 
enabling one to define a (p c )eff via extrapolation or just as the real part of the nearby 
endpoints. Our results have the property that these values of p c or (p c )eff he in the 
interval p G (1,2]. It is interesting to discuss the relation between this result and the 
Brown-Colbourn theorem [15] that R(G,p) does not have any real zeros in the region 
outside of the set U (1, 2]. In general, the property that a polynomial has no zeros in 
an interval of the real axis does not exclude the possibility that as the degree of this 
polynomial goes to infinity, a resultant continuous accumulation set of zeros B crosses 
the real axis in this interval. A well-known case where this phenomenon does occur is 
provided by the Potts model partition function. This is a polynomial in q and a = e K 
with non-negative coefficients. Thus, for fixed positive q, for either sign of the spin-spin 
coupling J, this partition function cannot have any zeros on the positive real a axis. 
However, for the infinite limit of the square lattice, the locus B in this case does cross 
the positive real a axis at the point a c = 1 + y/q, which serves as the phase boundary 
between the paramagnetic and ferromagnetic phases occupying the respective intervals 
1 < a < a c and a c < a < oo. 

For all of the lattice strip graphs for which we have calculated R(G,p) and B, we find 
that B is independent of the longitudinal boundary conditions, although it depends 
on the transverse boundary conditions and the type of lattice. This requires that 
the dominant a in the case of cyclic strips must arise from the subset of the cc's 
with coefficient that coincides with the set of en's for the corresponding strip with 
free longitudinal boundary conditions. The property that, for a given type of lattice 
strip, the locus B is independent of the longitudinal boundary conditions is quite 
different from what we found for the corresponding continuous accumulation sets of 
zeros of chromatic polynomials of lattice strips [30], [42] -[50], which do depend on both 
longitudinal and transverse boundary conditions. 

For all of the lattice strip graphs for which we have calculated R(G,p) and £>, we find 
that this locus consists of arcs (and in some cases a line segment) and does not enclose 
regions in the complex p plane. In some simple cases B is connected, but in general, it 
may consist of several disjoint components. Again, this is quite different from what we 
found for the loci B for chromatic polynomials of lattice strips; for those, a sufficient 
(not necessary) condition that this locus encloses areas is that the strip has periodic 
longitudinal boundary conditions [42, 46, 30]. 

In all of the cases for which we have calculated exact results, whenever B consists of 
an arc of a circle, then this circle is centered at p — 1. More generally, even if B is not 
an arc of a circle, it is roughly centered around p — 1. For families without multiple 
edges, we find that the component(s) of B is (are) roughly concave toward the left. 

The width L y — > oo limit of loci B for the infinite-length lattice strips is not necessarily 
the same as the continuous accumulation locus of the zeros of the reliability polynomial 
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for the two-dimensional thermodynamic limit L x — > oo, L y — > oo with L y /L x fixed to 
a finite nonzero constant. This is clear since for any finite L y regardless of how large, 
the infinite-length strip is a quasi-one- dimensional system, characterized by the ratio 
L y /L x = 0. Thus, for example, the Potts model on these strips has a paramagnetic- 
ferromagnetic (PM-FM) phase transition point only at K — oo, i.e., v = oo, while for 
2D lattices, it has such a PM-FM phase transition point at finite v. However, in our 
previous studies of continuous accumulation sets of loci of Potts model partition func- 
tions (e.g., [28]- [37] and our earlier work referred to therein on complex-temperature 
phase diagrams for Ising models) we have found that it is possible to gain some in- 
sight into certain features of the continuous accumulation set of zeros of the Potts 
model for the two-dimensional thermodynamic limit from studies of wider and wider 
infinite-length strips. 

In this context, we note that for the (infinite) square lattice, the PM-FM phase transi- 
tion point for the Potts model is given by v 2 = q [39] . This transition point is reflected 
in a nonanalyticity in the free energy of the Potts model at the corresponding tempera- 
ture variable v. As we have discussed and studied in earlier papers, e.g., [52, 53], when 
one generalizes v to complex values, one sees that a singular point for a real physical 
value of v, such as the PM-FM transition point, is associated with the fact that a 
complex-temperature phase boundary crosses the real physical v axis at this point. At 
the crossing points there are physical or complex-temperature singularities in the free 
energy [52]- [56]. The boundaries on B (here as a function of v for fixed q, but more 
generally, also as a function of q for fixed v) arise as the continuous accumulation set 
of the complex-temperature zeros of the partition function, analogous to the origin of 
the locus B studied here as continuous accumulation set of the zeros of the reliabil- 
ity polynomial. Connected with this, the physical paramagnetic, ferromagnetic, and 
antiferromagnetic (AFM) phases have complex-temperature extensions, and these are 
separated from each other by the phase boundaries B in the v plane. We recall that 
for the Potts model on the square lattice, the PM and FM phases occupy the intervals 
< v < y/q and y/q < v < oo, respectively. Analytically continuing q to real values, 
using eq. (3.6), we infer that in the limit q — > that yields the reliability polynomial, 
as given by eq. (3.11), the PM phase therefore contracts to a point at v — 0, which 
corresponds, via the relation p — v/(l + v) in eq. (3.12), to the point p — 0, while the 
FM phase expands to occupy the interval < v < oo, corresponding to the interval 
< p < 1. As is evident, e.g., in the square-lattice Ising model and more generally 
for the 2D g-state Potts model, the complex-temperature extension of the FM phase 
includes part of the negative real v axis extending to v — — oo (as well as the outerlying 
region \v\ — > oo in the v plane). This semi-infinite interval of large negative real values 
of v is mapped, via the transformation p — v/(l + v), to a finite interval of values of 
p > 1. Points in the outer region v = Re %e with R » 1 are mapped via this transfor- 
mation to points in the vicinity of p = 1, namely p ~ 1 — Rr l e~ %e . Since (in the case 
q — > relevant here) all of the points mentioned above in the v plane, namely (i) the 
interval < v < oo, (ii) the interval of large negative real v, and (iii) the semi-infinite 
annular region v = Re td with R — > oo are in the complex-temperature extension of 
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the FM phase of the Potts model, it follows that the images of these sets of points 
are all in the same analytically connected region of the p plane as defined by the locus 
B. We can draw on another source of information about this boundary as follows. 
Arguments have been given that for the Potts model on the (infinite) square lattice, 
the paramagnetic-antiferromagnetic (PM-AFM) transition occurs at a (physical) root 
of the equation v{y + 4) + q = [58]. Setting q = and v — p/(l — p) as in (3.11) 
and (3.12), one finds the solutions v = and v = —4, corresponding respectively to 
p — 0, the solution already found above, and p = 4/3. Hence for the square lattice, 
the image, under the mapping v — > p — v/(l + v), of the complex-temperature ex- 
tension of the FM phase includes the real interval < p < 4/3. The complementary 
intervals — oo < p < and 4/3 < p < oo would be part of the image under this 
map of an unphysical phase denoted the "O" (for "other" ) phase in Ref. [52] and our 
subsequent papers on complex-temperature phase diagrams. Using the fact that the 
complex-temperature phase boundary separates these various phases, we infer that for 
the reliability polynomial on the square lattice, in the thermodynamic limit, the zeros 
accumulate on a locus B that is a closed curve crossing the real axis at p = and 
p = 4/3 and separating the interior region from the exterior, which latter includes the 
intervals — oo < p < and 4/3 < p < oo. In an obvious nomenclature extending 
that of [52], we denote the regions in the interior and exterior of this closed curve as 
Rfm and Ro, indicating the correspondence with the complex-temperature phases of 
the Potts model. Just as the free energy is an analytic function in the interior of a 
complex-temperature phase, the function r that we have introduced in [14] and here in 
eq. (1.3) is an analytic function in the regions bounded by the locus B. One expects 
qualitatively similar results to hold for the regions of analyticity of r on (the thermo- 
dynamic limits of) lattice graphs of dimensionality d > 3. For the (infinite) square 
lattice, 

4 

Pc(sq) = - (7.1) 

We shall discuss how this compares with our exact results for p c values on infinite-length 
finite-width strips below. 

Similar reasoning can be applied to obtain inferences for the loci B and the correspond- 
ing region diagrams of analyticity of the function r on other (infinite) 2D lattices. For 
the triangular and honeycomb lattices, the PM-FM phase transition point is given by 
a (physical) root of the equations v 2 (v + 3) — q = and v 3 — 3qv — q 2 = 0, respectively 
[60]. Evaluating these for q = and v — p/(l—p) yields for the triangular lattice, the 
solutions p = and p = 3/2, so that 

p c (tri) = \ (7.2) 

For the honeycomb lattice, the same evaluation yields just the solution p = and 
hence does not determine p c (hc) (By the reasoning given above, one can infer that 
p c {hc) > 1). The fact that the point p = on B is common to all of the three lattices 
- square, triangular, and honeycomb - follows because it is the image under the map 
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(3.12) of the point v — K — 0. This point corresponds to infinite temperature, where, 
as one can see from eqs. (3.3) and (3.4), the spin-spin interaction does not contribute 
to Z(G,q,v), which reduces to for all of these lattices (indeed for any graph G). 

• For a given type of lattice strip (square, triangular, honeycomb, etc.) and a given 
set of transverse boundary conditions, we may consider the sequence of loci B for each 
width, L y , and inquire whether these approach a limiting locus as L y — > oo. Our results 
are summarized in Table 3 and are consistent with the hypothesis that such a limiting 
locus exists. As part of this analysis, we have investigated, for each type of strip family, 
how p c or ip c )eff depends on L y . We find that, for a given type of lattice strip and 
transverse boundary conditions, if a value of p c exists, i.e., if B crosses (intersects) the 
real axis, then p c decreases as L y increases. However, if, for a given type of lattice 
strip and transverse boundary conditions, one considers the set of values of both p c 
and (p c )eff (the latter for widths where B comes close to but does not cross the real 
axis), then the dependence on L y does not appear to be monotonic. For example, for 
the infinite-length limit of the strip of the square lattice with free transverse boundary 
conditions (and any longitudinal boundary conditions), for the width L y = 2, we have 
shown that p c = 4/3. Interestingly, this is the same value as inferred above for the 
infinite square lattice defined via the usual thermodynamic limit. For the square-lattice 
strip with width L y = 3, the locus B does not cross the real axis, but, as is evident 
in Fig. 2, the inner endpoints of the arcs on B are sufficiently close to the real axis 
that one can infer a (p c )eff, and it is (p c )eff — 1-335, which is slightly greater than 
the value of p c for the L y = 2 strip and the inferred value of p c {sq)- As is evident 
in Table 3, the values of p c are not as close to p c {sq) for the square-lattice strips of 
the corresponding widths with periodic transverse boundary conditions or self-dual 
boundary conditions, as compared with free transverse boundary conditions, although 
they are again consistent with approaching p c (sq) = 4/3 as the width increases. For 
the L y = 2 strip of the triangular lattice, we find the exact result p c = 3/2. Again, 
interestingly, this is equal to the inferred value of p c (tri) for the infinite triangular 
lattice defined via the thermodynamic limit. We have encountered this sort of situation 
before. For example, for the chromatic polynomial P(G,q) = Z(G,q,v = —1) and its 
asymptotic limiting function as |V| — > oo, W({G},q) (the degeneracy per vertex for 
the Potts model), we found that for the infinite- length strip of the triangular lattice 
with width L y = 3 and toroidal or Klein-bottle boundary conditions, q c (defined as the 
maximal point where the continuous accumulation set of zeros of P(G, q) intersects 
the real q axis) is q c = 4 [49], the same value as for the infinite triangular lattice. 
Similarly, for infinite-length limits of self-dual strips of the square lattice, we found 
q c = 3 [36, 37], the same value as for the infinite square lattice. 

As regards the behavior of B near p — 0, our results are also consistent with the 
possibility that for each type of lattice, as the strip width L y — > oo, the respective 
loci B have, on the left, complex-conjugate arc endpoints that come together and 
join at p — 0, as motivated by the discussion above for the infinite 2D lattices. The 
observed feature that this occurs only as a limit, in contrast to the crossing on the 
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right at the respective values of p c , can be understood via the relation (3.11) and the 
aforementioned fact that the infinite-length finite- width strips are quasi- ID systems as 
far as their thermodynamic behavior is concerned. Since, as we are noted above, the 
interior of B for the reliability polynomial on the 2D lattice corresponds to the image 
of the complex-temperature extension of the FM phase of the 2D Potts model (for 
q — > 0), and since the argument that this must be separated from other phases relies 
on the existence of nonzero ferromagnetic long range order (magnetization), it follows 
that for quasi-lD systems, where for a short-range interaction like that in (3.4), the 
standard Peierls argument shows that there is no long range order, there is no FM 
phase. This yields a plausible explanation of the features that we observe in our exact 
solutions for the infinite-length, finite- width strip graphs, that the loci B do not enclose 
regions in the p plane. 



r 




Figure 17: Comparison of asymptotic reliability per vertex r for the L x — > oo limit of several lattice strips. 
From bottom to top, the curves refer to (a) sq[l x oo, FBC y ) with d e ff — 2, (b) hc[2 x oo, FBC y ] with 
d e ff — 2.5, (c) sq[2 x oo, FBC y ] with d e ff — 3, (d) sq[2 x oo, PBC y ] with d e ff — 4, (e) Gd[1 x oo], same r 
as for tri[2 x oo], with d e ff = 4, (f) Gd[2 X oo], with d e ff = 4. 



8 Conclusions 

In this paper we have presented exact calculations of reliability polynomials R(G,p) for 
lattice strips G of fixed width and arbitrarily great length with various boundary conditions. 
We have introduced the notion of a reliability per vertex, r({G},p) and have calculated this 
exactly for the infinite-length limits of various lattice strip graphs. We have also studied the 
zeros of R(G,p) in the complex p plane and determine exactly the asymptotic accumulation 



54 



Table 3: Properties of reliability polynomials and resultant loci B including p c or {p c )eff f° r various lattice 
strip graphs. The notation sqdl and D refer to the self-dual strips of the square lattice. 
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set of these zeros B, across which r({G}) is nonanalytic. We have observed and discussed 
several general features of the r functions and the loci B for these families of graphs. 
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